Structural anisotropy and orientation-induced Casimir repulsion in fluids 
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In this work we theoretically consider the Casimir force between two periodic arrays of nanowires 
(both in vacuum, and on a substrate separated by a fluid) at separations comparable to the period. 
Specifically, we compute the dependence of the exact Casimir force between the arrays under both 
lateral translations and rotations. Although typically the force between such structures is well- 
characterized by the Proximity Force Approximation (PFA), we find that in the present case the 
microstructure modulates the force in a way qualitatively inconsistent with PFA. We find instead 
that effective-medium theory, in which the slabs are treated as homogeneous, anisotropic dielectrics, 
gives a surprisingly accurate picture of the force, down to separations of half the period. This 
includes a situation for identical, fluid-separated slabs in which the exact force changes sign with 
the orientation of the wire arrays, whereas PFA predicts attraction. We discuss the possibility of 
detecting these effects in experiments, concluding that this effect is strong enough to make detection 
possible in the near future. 



I. INTRODUCTION 

Casimir forces are usually attractive interactions mea- 
surable at small separations, but recent theoretical 
works [THE] have predicted a variety of situations in 
which these forces can be modified by using complex mi- 
crostructures. In addition, by utilizing different choices 
of materials the Casimir force can be changed in both 
magnitude [9 and sign [T0HT2] . However, with some ex- 
ceptions [7 , 8j, the qualitative aspects of these effects can 
be explained through a combination or competition of 
forces calculated using some form of the proximity force 
approximation (PFA, a common heuristic description of 
the Casimir force as pairwise interactions between par- 
allel surface patches) [13 j. It is therefore of interest to 
consider situations in which geometry not only allows 
the Casimir force to be modulated (e.g., by reducing or 
changing the sign of the force), but also creates effects 
that cannot be accounted for by PFA. 

In this article, we introduce and examine a geometry 
that exhibits both of these qualities. We examine the 
configuration shown in Fig. [I] consisting of two identical 
microstructured slabs consisting of periodic arrays of di- 
electric nanowires. We compute the exact Casimir force 
for this geometry using a combination of existing scatter- 
ing theory techniques [TIHTS] , to be described below. We 
find that the microstructure of the slabs leads to a num- 
ber of interesting qualitative effects: the force between 
the slabs can be dramatically modulated by rotating the 
two slabs at fixed surface-surface separation. In vacuum, 
the force between the slabs can be significantly reduced 
(more than halved as the rotation angle changes), and if 
the slabs are immersed in a fluid, the sign of the force 
can flip. It turns out that, even at moderate separations 



d/a ~ 1 5 PFA (an uncontrolled approximation) cannot 
capture either of these effects. Specifically, consider the 
behavior of the Casimir force between the slabs as they 
undergo lateral translations y (at y = the slabs are 
mirror-symmetric) and rotations 6. At = 0, the slab 
microstructures are aligned as in Fig. [I] (upper), while 
for — it/ 2 the slabs are crossed, as in Fig. [I] (lower). 
As y and are varied, the slab surface-surface separa- 
tion d is kept fixed. From a simple geometric argument, 
it is clear that for vacuum-separated wires PFA predicts 
the following bound on the Casimir force between two 
identical slabs: 

^aligned, y —a/ 2 ^ -^crossed ^ -^aligned, y=0 (1) 

(a positive force is attractive). Here Crossed denotes the 
force when 6 = n/2. This bound is insensitive to the de- 
tails of the exact PFA used, and simply relies on the fact 
that the wire surface-surface separations are minimized 
at 6 = 0,y = and maximized when = 0, y = a/2. 
Although a PFA prediction must be valid as d/a — > 0, 
for the systems examined in this work we find that the 
bound Eq. ([!]) is violated even at moderate distances 
d/a ~ 1. For vacuum- separated metallic/dielectric (e.g., 
gold or silica) nanowires and no substrate, this implies 
that changing the orientation of the wire arrays (i.e., 
their geometry) plays a stronger role in reducing the 
force than simply reducing the pairwise surface-surface 
separation between the slabs. In another case to be dis- 
cussed, in which the wires are gold, the substrate silica, 
and the fluid ethanol, the bound Eq. ([!]) is also valid. 
However, we find from exact calculation that while PFA 
predicts F y=a / 2 < and Crossed > 0, the opposite is 
in fact true, i.e., aligned slabs are always attracted to 
each other, while crossed slabs are repelled. Therefore, in 
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FIG. 1: Configuration examined in the text. Each slab con- 
sists of a single periodic array of cylindrical wires embedded in 
a semi-infinite substrate; the goal is to create for each slab an 
effective anisotropic medium. Of interest here is the change in 
the force from when the wires are aligned (Top) and when they 
are crossed (Bottom). When the medium between the plates 
is a properly chosen fluid, we show that the force switches 
sign as a function of orientation for certain separations. 



our case at d/a ~ 1 the sign of the force can therefore 
be modulated by rotating the slabs (but not translating 
them), an interesting possibility for experiments that we 
discuss later. 

To understand these effects qualitatively, we find 
that a much more suitable framework is the effective 
medium approximation (EM A), in which the microstruc- 
tured slabs of Fig. [I] are describable as an homoge- 
neous medium with a given anisotropic permittivity ten- 
sor s= diag (s\\,s±,Sj_) (here and below, e, e\\ and e± 
will be used to denote effective, homogenized permittiv- 
ities). We are able to compute the effective e tensor for 
this configuration from our scattering method [17 , and 
we find that for the gold-silica-ethanol configuration, the 
deduced effective-medium parameters follow an ascend- 
ing sequence e± < enuid < This ascending sequence 
is known to lead to repulsion [18] for uniform, isotropic 
materials, and we find below that the effect is also present 
for anisotropic materials and individual polarizations (al- 
though in this case, as discussed below, an ascending se- 
quence is not sufficient for repulsion). One expects the 
EM A to be strictly valid as d —> oo, however in our case 
it turns out, somewhat surprisingly, to be qualitatively 
accurate down to d/a < 1 (a mathematical proof of why 
the EM A holds at such short separations, at least for 



crossed slabs, is given in Appendix [B| . Further, we find 
that for realistic materials (gold wires, silica substrate, 
immersed in an ethanol fluid) and geometry parameters, 
this orientation induced repulsion at fixed surface-surface 
separation between the two slabs holds down to separa- 
tions d/a ~ 0.5. 

Previous works [TTJ [19] have discussed the possibility 
of using geometry to create artificial constitutive rela- 

tions (e.g., effective permeability | M | > | s | [under 
a suitable matrix norm] to obtain Casimir repulsion or 
chirality n ^ to achieve force reduction). However, it 
has been rigorously proven [20] that no such media with 

M> £ can be constructed from metallic/dielectric con- 
stituents so as to exhibit repulsion in vacuum, and com- 
putations involving the exact microstructures [21] have 
shown that chirality effects are only present at separa- 
tions so large that they cannot conceivably be detected. 
Therefore, while in some circumstances EMA can be a 
useful qualitative guide (and rigorously accurate in cer- 
tain limits), it must be used with caution — ideally, as 
a supplement to exact calculations. Orientation depen- 
dence (and the resulting Casimir torque) between slabs 
has previously been considered for two birefringent plates 
with weak anisotropy [22H24] , and for corrugated metal- 
lic plates [2]. Although the torques in these systems are 
in principle measurable, the change in the force with ori- 
entation is small, and the forces are always attractive. 
Ref. [25] showed a large orientation dependence for theo- 
retical uniaxial conductors, and suggested a possible re- 
alization via nanowire arrays, but without calculation in 
the latter case and without changing the sign of the force. 

In the present work, we analyze this effect 
for periodically-patterned, vacuum-separated suspended 
membranes (which form a potentially promising medium 
for Casimir force measurements [26]) and for gold wires 
on silica substrate, immersed in ethanol. In both cases, 
we find Eq. (TJ is violated down to d/a ~ 1. In the former 
case, we find that the force is 70% lower for crossed slabs 
compared to aligned slabs (for any y) at large separations, 
while at separations comparable to the unit cell (e.g., 
100 nm) there is a more modest, but still significant, 30% 
reduction. For the latter system, we find that aligned 
slabs are attractive for all y down to d/a ~ 0.5 while 
crossed slabs are repulsive in this range. These examples 
demonstrate a system in which effective-medium theory 
is correct and gives predictions that differ from PFA in 
a highly nontrivial way. Given that nanowire arrays be- 
low 15 nm can be fabricated with current technology [27] 
and that for these dimensions the predicted effects occur 
down to sub-lOOnm length-scales, these effects should be 
experimentally accessible. After presenting these results, 
we argue how they may be experimentally detectable 
(assuming suitable fabrication techniques) and estimate 
force magnitudes in the hypothetical case where one of 
the slabs is replaced by either a sphere or cylinder with 
a wire pattern stamped on its surface. 
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II. METHOD 

In this work, we perform Casimir force calculations 
with a semi-analytic scattering method using a combina- 
tion of results from Ref. [T4lfT6l which efficiently computes 
the exact Casimir force between periodic arrays with one 
dimensional of translation-invariance. Our implementa- 
tion differs somewhat from previous scattering methods 
(e.g., Refs. 28, 29) in that it is particularly well-suited to 
unit cells with objects of compact cross-section, such as 
the wires in Fig. [I] (rather than, e.g., extended rectangu- 
lar gratings). While we are primarily interested in objects 
of circular cross-section in the current work, the present 
method can be extended in a straightforward manner to 
treat unit cells with arbitrarily shaped compact objects 
using an existing boundary-element method [30j [31] . In 
addition, we have checked the results with a brute-force 
finite-difference time-domain (FDTD) [32j [33] method 
and found good agreement, and additionally these com- 
putations show similar results for square wires. Also, the 
effects are not significantly different at zero temperature 
and 300 K, so we work exclusively in the former limit. 

The two plates are separated from each other by a 
distance d in the x direction, and the slab is termed y- 
directed if the wires are along the y axis (similarly for z) . 
The slabs are aligned if both slabs are y- or z-directed, 
and crossed if they are not. The zero-temperature 
Casimir interaction energy for a ^-directed slab displaced 
by x from a q' directed slab (<?, q' = z) is: 
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(2) 

where Ri 5 2 are the scattering matrices of planewaves for 
the two slabs, U(x) is the planewave translation matrix 
for relative displacement x between the slabs, and the in- 
tegral of transverse wavevector components k yi k z is over 
the first Brillouin zone. See Ref. [l6|for a detailed deriva- 
tion and a partial review of precursors [34-36 . For the 
present work, we require an efficient method of comput- 
ing the scattering matrices from periodic arrays, which 
we describe in more detail in Appendix | A 1[ 

For sufficiently large distances, the relevant frequencies 
and wave- vectors are unable to probe the structure of the 
arrays and consequently an effective medium approxima- 
tion should produce good results. In our case, such an 
effective medium should have a much larger conductivity 
in the direction of the wires as compared to the other 
orthogonal directions even in the static limit (this being 
one of the motivations for the current configuration), so 
it is clear that an anisotropic EMA is called for. Fortu- 
nately, Casimir interactions between anisotropic homoge- 
neous media have been studied by several authors [T71 I37]- 
[40], allowing us to build upon their resul ts. W e use the 
method of Ref. [171 (outlined in Appendix A 2 ) to obtain 
the scattering matrices assuming a known permittivity 
e . In the subsequent analysis, we also require the inverse 
procedure: given the scattering matrices (computed us- 



retrieve the best-fit EMA e . This procedure quickly be- 
comes intractable if arbitrary e and M are allowed. In- 
stead, we assume /x^ = Sij and s xx =e yy exhibit no de- 
pendence on aside from its polarization. For —> 
the scattering matrices for each polarization reduce to the 
standard Fresnel formula for reflection off isotropic inter- 
faces. Their inversion then yields the effective dielectric 
tensor: 



6\\(i£) = Sfluid 



1-R 



1 + R 



e±(i£) = £fluid 



l + R± 
1-R± 



(3) 



III. RESULTS 



In the configuration of Fig. [T] the wires have radius 
r = 0.3 a and period a; we take the wire centers to be 
in the substrate and the wire surface to be tangent to 
the substrate surface to maximize the wire-wire interac- 
tions. The permittivities of wires, substrate, and fluid 
are respectively £2^0 and €3(1%) for imaginary 

frequency £. The materials used for the wires, substrate, 
and fluid are gold, ethanol, and silica, respectively. For 
gold, we use a plasma model with £ p = 1.36 x 10 16 rad/sec 
(the addition of a small loss term does not change the 
results significantly). For silica and ethanol we use stan- 
dard oscillator models [4lJ [42] . 



A. Vacuum- separated slabs 

We first compute the Casimir forces when the inter- 
vening medium is vacuum (e 3 = 1), comparing the force 
^aligned (v) f° r aligned slabs with the force F crosse d for 
crossed slabs. In the aligned case, the force can depend 
on y (leading to a lateral component of the force), and 
in this case ^aligned refers only to the normal compo- 
nent of this force. In the crossed case = ir/2 there 
is no ^/-dependence. Although there is no sign change 
for vacuum-separated slabs, this configuration is of in- 
terest because it can be fabricated as a single suspended- 
membrane structure for each orientation [26] , and may be 
easier to work with than a fluid system. Figure [2] shows 
results for wires composed of perfect metal, gold, and 
heavily doped silicon, the latter being more conventional 
for fabrication. We plot the ratio F crosse d / ^aligned (y) (the 
shaded regions indicate the full range of this ratio as 
y is varied), which serves two purposes. First, it indi- 
cates the required relative accuracy in a force measure- 
ment needed to discern the orientation-dependence of the 
force in an experiment. Second, it indicates the tran- 
sition from the PFA to the EMA regimes via Eq. ([!]): 
if the force is determined as a sum of pairwise interac- 
tions, then the total force is maximized when the pair- 
wise surface-surface distances are minimized. A simple 
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geometric argument shows that the net distance is max- 
imized for = 0,y = a/2. On the other hand, this dis- 
tance is minimized for — 0, y = 0. The bound Eq. ([I]) 
then follows. It turns out that the force is always at- 
tractive in this situation, so we have the further bound 
^aii g ned(^ = a/2) > 0. Therefore, PFA predicts that 
there is a range of y such that F crossed (y) /^aligned > 1- By 
contrast, EM A predicts the inequality F crosse d < ^aligned 
(there is no ^/-dependence in this approximation). This 
inequality stems from the following scattering-theory ar- 
gument: for d/a ^> 1, the exponential suppression of 
k y , k z in U [16] implies that the force is domi- 
nated by the scattering of planewaves at normal incidence 
(k y = k z = 0). For normal incidence, the reflection ma- 
trix R(*) is anisotropic but diagonal in polarization, and 
can be computed from an effective anisotropic dielectric 
tensor: e (it;): 

*e (ig) = diag (e xx , e yy , e zz ) K « diag (R\ \,R±) (4) 

where R\\ and R± are the matrix elements for incident 
light polarized parallel and perpendicular, respectively, 
to the wires. For normal incidence, e xx does not con- 
tribute, and we define e\\ = e zz and e± = e yy . Due 
to the high permittivity of the wires and the low per- 
mittivity of the substrate, we expect \R±\ <C \R\\\ and 
R±R\\ > [note that R(i£) is real]. When the wires are 
aligned, the relevant product of scattering matrices in the 

energy integrand is (^R^ + R 2 ^ ~ R^ , and when they are 

crossed the term is 2R\\R ± , implying F crosse d < ^aligned- 
From Fig. |2j this bound is clearly violated in doped Si 
for d/a < 0.6, and for other materials at smaller d. Note 
that this includes a region d/a < 2 for which there is 
still a strong ^-dependence of this ratio, indicating that 
while EMA is not strictly valid, it still has more predic- 
tive power than PFA. 

As d/a — >• oo, the force ratios should approach a con- 
stant determined by e (£ — »• 0): this ratio is ~ 0.29 for 
all three materials. However, this limit is approached 
very slowly (0(<i -1 / 2 ) for perfect metal/gold wires). The 
EMA ^-dependence is due both to constituent mate- 
rial dispersion and an effective geometric dispersion in 
e (i£) induced by the geometry of the wires. As £ — > 0, 
£|| —> e±(i£) and e± —> €2(1^)- The former limit follows 
from the l/£ 2 divergence of S\ as £ — >• (if we add a 
dissipative term to £1, the divergence is only l/£ but 
this crossover occurs at frequencies too low to affect the 
present results). The latter limit follows from the fact 
that the static polarizability of the wires in the transverse 
direction is finite, implying that their contribution to the 
scattering amplitude vanishes as £ — » 0. Therefore, as 
i — » 0, £±_ depends only on the substrate, which has semi- 
infinite extent; we have checked this relation numerically 
and found good agreement, lim^oo Crossed /^aligned can 
then be obtained by computing the force between two 
anisotropic plates, with s\\ = 00 and e± = £2(0). 

For all d/a ^> 1.5, the orientation-dependence of the 
force is quite strong (^> 30%). In this range, the absolute 
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FIG. 2: Ratio of the Casimir force between crossed wire arrays 
(^crossed) over aligned wires (i 7 kiigned(2/)) between vacuum- 
separated nanowires of perfect metal, gold, and doped silicon 
(dopant density 10 20 /cm 3 ). a — >• 00 corresponds to using non- 
dispersive e(i£ — ► 0) for all materials (because £ is in units of 
c/a). Shaded regions indicate the range of values that ^aligned 
takes over all relative lateral displacements y of the wire cen- 
ters, with their higher (lower) boundary indicating the force 
for y = (y = a/2). PFA predicts, via Eq. |l]), that the max- 
imum of this ratio should exceed 1, which is only true for Si 
at d/a ~ 0.6, and for gold/perfect metal at lower separations. 
Inset: details of the configuration examined. 



pressure is approximately 2% of the corresponding pres- 
sure between two homogeneous perfect metal plates. As 
differences of this magnitude between vacuum- separated 
plates have been measured on the 100 — nm length scale, 
we are hopeful that this departure from PFA can be de- 
tected experimentally. 



B. Fluid-separated slabs and tunable repulsion 

In the previous example the force is always attractive, 
because £i(i£), £2(^0 > e 3^0 f° r a U £ ED]- However, 
if the medium between the slabs is a fluid such that 
£i(*0 > £3(*0 > £2(*0> the above EMA analysis pre- 
dicts that the force in the crossed configuration should 
be repulsive as d/a — » 00 since by the Fresnel formula 
at £ = 0, R±R\\ < 0, leading to the well-known Casimir 
repulsion for an ascending sequence of e for each polar- 
ization [TT]. By contrast, in the aligned configuration 
the force is attractive when y = due to mirror symme- 
try [43], and the intuition that this attraction holds for 
non-zero y is confirmed by calculation. However, strong 
corrections to this argument occur at finite separations: 
waves with nonzero k y and k z contribute, and are not in 
general polarized along the y or z axes. These waves will 
couple to both e\\ and e± on each reflection. The sign 
of the resultant reflection coefficient will usually be the 
same for both z- and ^-directed slabs, leading to attrac- 
tive contributions. Further, the microstructure can also 



5 



lead to significant corrections, possibly eliminating the 
effect for separations comparable to the unit cell a [2T] . 
As such separations are necessarily required for experi- 
ments, we require exact results to verify that this effect 
persists in experimentally accessible regions. 

The exact results are shown in Fig. [3] as in Fig. |2j 
shaded regions show the ^/-dependence. The results show 
a clear attractive-repulsive transition as the slab orienta- 
tion is varied between the aligned and crossed configura- 
tions. This effect persists for both the ideal case of perfect 
metals (a — >• oo) and dispersive materials at a = 100 nm; 
both show qualitatively similar behavior. For crossed 
slabs, the repulsion is fairly flat over a ~ 40 nm range. 
At first sight, it is temping to ascribe the repulsive force 
observed for crossed wires to PFA-line interactions be- 
tween opposing areas where metallic wires face dielectric 
substrate and thus feel a repulsive force. However, if this 
were the case, then by the same argument (and Eq. ([T])), 
^aligned (y = a/2) should exhibit an even stronger repul- 
sion. This is clearly not the case in Fig. |3| for d/a > 0.5; 
rather, the results for repulsion are consistent with the 
EM A argument given above. To estimate the magni- 
tude of this repulsion, for a = 100 nm the repulsive pres- 
sure is approximately 1% of the pressure between paral- 
lel perfectly-conducting plates separated by a compara- 
ble distance d in vacuum, e.g, ~ 0.13pN//im 2 . Repulsive 
forces in fluids on this order of magnitude have been mea- 
sured [IT], therefore these forces are potentially within 
reach of current or near-future measurement techniques 
(the force in more realistic sphere-plate and c ylinder- 
plate geometries will be considered in Sec. HID below). 
Interestingly, when the wires are crossed and d varies 
there is an attractive-repulsive transition at a critical 
separation d — d c . This leads to an unstable equilib- 
rium with respect to d. (The crossed configuration is 
always unstable with respect to orientation.) This tran- 
sition, as mentioned above, is attributed to the eventual 
dominance of attractive forces as separation goes to zero. 
However, we will see below that such a transition is pre- 
dicted by EMA as well, and is therefore not due entirely 
to proximity effects. 

Before continuing, it is interesting to see what the PFA 
prediction for the force in this system is. We simplify 
matters by computing the PFA assuming instead wires 
of square r x r cross-section, with centers a distance r/2 
beneath the substrate surface (by minimizing the amount 
of curved surface, we expect this to maximize the range 
of separations for which PFA is accurate). Using FDTD 
computations [32 , 33 , we have found that this system ex- 
hibits behavior qualitatively similar to that of Fig. [3] The 
results for a = 100 nm are shown in Fig. [4] (taking a — »• oo 
yields a qualitatively similar curve). Inspection reveals 
that the PFA prediction exhibits the opposite behavior 
for the sign of the force as the full numerical calcula- 
tions: the PFA force is exclusively attractive for crossed 
wires, but for aligned wires shifted by y = a/2 the force 
becomes repulsive. 

Although PFA fails qualitatively and quantitatively to 
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FIG. 3: Casimir pressure for the same parameters as the inset 
to Fig. [2] with a silica substrate and ethanol between the two 
slabs. Positive values indicate an attractive force, negative 
values (shaded) repulsive. Both a — » 00 and a = 100 nm are 
shown. As in Fig. [2] shaded regions denote the range of y- 
displacements. When the wires are aligned (blue), the force 
is always attractive, but when the wires are crossed (red), the 
force turns repulsive at a critical separation d > d c , which 
depends on a. 
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FIG. 4: Forces computed for the configuration of Fig. fusing 
the Proximity Force Approximation (PFA), for a = 100 nm. 
For convenience in applying the PFA, square wires of width 
r are used instead of circular wires (the exact results are not 
qualitatively changed). We see that the behavior of Fig. [3] 
is not qualitatively captured: the force is always attractive 
for crossed wires, and the bound Eq. |l]) is strictly satisfied. 
Similar results are obtained for a — » 00. 



predict the results of Fig. [2] and Fig. [3j we have not 
yet examined the quantitative accuracy of EMA. In the 
next section, we will examine a description of the force 
in terms of the EMA, rigorously valid in the regime 
(d/a — » 00). We will see that EMA gives qualitatively 
correct predictions for the magnitude and sign of the 
forces down to separations comparable to the unit cell 
size. 



6 




FIG. 5: Effective-medium parameters for the system of Fig. [3] 
derived from scattering data for a length scale a = 100 nm. 
£\\ is the effective dielectric parallel to the wires (black), and 
e± the effective dielectric perpendicular to the wires (blue). 
The dielectrics of the constituent materials gold, silica, and 
ethanol are shown for reference. These dispersions will be 
used to compute the Casimir force in the effective-medium 
approximation below. 



C. Comparison with EMA 



FIG. 6: Results for the force computed in the Effective- 
medium Approximation (EMA) using an effective homoge- 
neous, anisotropic permittivity tensor ?, via the method de- 
scribed in Appendix |A 2| As in Fig. [3] results for both a 
infinite (i.e., non-dispersive materials) and a = 100 nm are 
shown, for both crossed and aligned configurations. Although 
quantitatively inaccurate over the range plotted, the EMA 
does qualitatively capture the behavior of the force a opposed 
to PFA. Also shown are d c , the location of the attractive- 
repulsive transitions for the crossed wires (d c = 0.18a for 
a — >• oo and is not shown). 



In this section, we examine the extent to which the 
EMA analysis given above in Sec. |III A| predicts the cor- 
rect results, using the method of Appendix |A 2| We use 

a simplified EMA, assuming that M = 1 and e (i£) 
is isotropic in the plane perpendicular to the wires. In 
this case, e\\ and e± are the only EMA parameters in- 
volved; these are retrieved by applying Eq. (|3| to the 
scattering matrices computed from the exact structure. 
The EMA-predicted force is then computed from the 
method of Appendix |A2| A plot of the retrieved pa- 
rameters for a = lOOnm is shown in Fig. [5] Shown for 
reference in Fig. [5] are e(i£) of the constituent materials 
ethanol, silica, and gold. As expected, e± approaches 
^silica f° r l° w while for higher £ the wires increase the 
effective permittivity. In e\\, we see a significant geomet- 
ric dispersion: in fact, at low £, ~ 1 + (£p/£) 2 ? 
where ^ p ~ 0.36 £ p is a new, effective plasma frequency 
for the gold wires. This is similar to an argument pre- 
sented in Ref. |44j where the effective dielectrics of square 
arrays of wires along the wire axis can be described by 
a plasma model with a reduced plasma frequency (as we 
have only included a single row of wires, Ref. H4js result 
cannot be directly applied, but its basic idea remains). 
A similar effect holds when the wires are perfect conduc- 
tors, where = 7.9. This accounts for the ^-dependence 
in the EMA regime of Fig. [2] the geometry introduces 
an effective length scale into the system, given by £ p . 

We use the retrieved £|| and e± to compute R and 
hence the Casimir force for the fluid-separated geometry. 
Computation of the force with the EMA parameters for 



the vacuum-separated slabs in the limit d — » oo gives 
agreement with the results of Fig. [2] The results for the 
fluid-separated case are shown in Fig. [6] We find that, as 
opposed to PFA, EMA gives qualitatively accurate (i.e., 
the same order of magnitude) predictions for both the 
magnitude and sign of the force. In particular, EMA 
predicts attractive-repulsive transitions at some d = d c , 
where for d < d c the force for crossed slabs is attractive. 
The fact that is a rapidly decreasing function of £, 

while e± increases, indicates that within the EMA the 
force should receive attractive contributions from higher 
£ (which dominate at small separations). Further, the 
decrease of both e\\ and with decreasing a also con- 
tributes to a reduction in repulsion. This explains why 
the repulsion for a = 100 nm is lower than for a —> oo. 
With our simplified EMA the predicted values of d c are 
very inaccurate: d c ~ 0.18 a for a = oo and d c ~ 1.05 a 
when a = 100 nm. However, the ^-independent regimes 
of Figs. 2-3 suggests that some EMA must be valid in 
those regimes, and in this case a more accurate EMA 
would involve both e xx ^ e± as well as an anisotropic 

effective A more general model of dispersion would 
allow for a translation-invariant but k-dependent permit- 
tivity. Fitting an effective s (k) is somewhat compli- 
cated in this framework; rather, to explore the validity 
of this "specular" approximation we compute the force 
with all non-specular (i.e., terms not conserving k^+G_L, 
where Gi is a reciprocal lattice vector) terms in the scat- 
tering matrix R removed. We find the surprising result 



that, while this specular approximation does not give sig- 
nificantly different predictions for the force for aligned 
wires, it gives much more accurate predictions for d c and 
the magnitudes of the force for crossed wires. This indi- 
cating that non-specular scattering events are suppressed 
when the wires are crossed, even at relatively small sepa- 
rations d/a ~ 0.5. We give a rigorous proof of this result 
in Appendix [Bj using a recently- developed diagrammatic 
expansion for the Casimir energy [45] . 



D. Forces for other geometries 

Parallel-plate configurations involving suspended 
membrane structures show potential for new sets of 
experimental Casimir force measurements for vacuum- 
separated geometries [26] ; such measurements detect the 
Casimir force (or force gradient) through a shift in either 
the optical spectrum or the resonance frequency of the 
upper membrane. However, to measure the sign change 
in the force predicted here a measurement between 
fluid-separated objects must be performed, where fluid 
damping is prohibitive. Instead, a force measurement 
involving an object (e.g., a sphere) mounted on an 
atomic force microscope (AFM) tip pj] is more realistic. 
In this case, alignment issues favor the use of one 
spherical and one planar object [46JH7], rather than the 
two planes considered here. For our purpose, a pattern 
of wires similar to Fig. [I] would be need to be stamped 
on the surface of the object. In this section, we give 
force predictions both sphere-plate and cylinder-plate 
geometries. The latter case has more difficult alignment 
issues, but still simpler than plate-plate alignment, and 
in this case the Casimir force is much larger than for 
sphere-plate setups [48]. For both cases, the radius of 
curvature R is many orders of magnitude larger than 
the surface-surface separation d (e.g., R = 200 /im and 
d ~ 100 nm). Because d, a <C R, it is appropriate to 
use a hybrid PFA/exact method in which each unit of 
surface on the sphere/c ylinder feels the exact pressure 
(as computed in Sec. |IIIB ) between a plate-plate 
configuration of the same surface-surface separation. 
The result is asymptotically exact in the limit R — >• oo 
for fixed d and a. 

Results are shown in Fig. [7] for R = 200 /im, and are 



plotted in units of pico-Newtons. The sphere-plate force 
has a peak repulsion of ~ 10 pN, which although below 
the detection limit of current measurement techniques in 
fluids, may be observable in the near future. We also note 
that, although the presence of Casimir repulsion may be 
experimentally challenging, Eq. ([I]) can be verified with 
much less sensitivity (e.g., ~ 50 pN at d = 50 nm). For 
the cylinder-plate, we take a length L = 16 jam in or- 
der for the force to be comparable in magnitude to the 
sphere-plate force. However, it is clear that if we instead 
take L ~ R, we can obtain both repulsive and attractive 
forces on the order of 10 3 pN, well within current exper- 
imental detection ranges. Therefore, this system seems 




Wires crossed 



50 100 150 

Surface-surface separation (nm) 

FIG. 7: Forces for sphere-plate (red) and cylinder-plate con- 
figurations, for a — 100 nm. In both cases, the wire pattern 
of Fig. [I] is stamped onto the surface of both objects, which 
have radii of curvature R = 200 /im. The cylinder length is 
chosen to be L = 16 /im so that the two force curves have 
comparable magnitude; however the cylinder-plate force is 
proportional to L and so will be much larger in magnitude 
for a more realistic L ~ R. 



most attractive for detection of an attractive-repulsive 
transition with orientation in a fluid if techniques similar 
to those of Ref. [49] to align the long axis of the cylinder 
relative to the plane can be extended to fluid-separated 
objects. 



IV. CONCLUSIONS 

We have presented an example system, consisting of 
microstructured slabs, in which PFA fails at relatively 
small separations. Instead, EMA qualitatively describes 
the behavior of the Casimir force, including the case in 
which the force can be switched from attractive to re- 
pulsive as the slabs are rotated. We have discussed the 
prospects for detecting these forces in experiments. One 
issue arising in an experiment, not discussed here, is 
the increased complexity of the electrostatic calibrations. 
Conceivably, these complications can be eliminated in 
the fluid-separated case by the addition of electrolytes to 
the fluid. This should also reduce or eliminate the effect 
of accumulated surface charges on the gold-silica inter- 
face, due to contact potentials [50j[5T. For the case of 
vacuum-separated doped silicon, a more involved calibra- 
tion procedure is required. 

Although not examined here, we also note that orient a- 
tional attractive-repulsive transitions may also occur us- 
ing naturally anisotropic materials. We have confirmed 
this for the case of lithium niobate slabs immersed in 
ethanol at zero temperature; although here d c > 7/im 
(implying that finite-temperature effects must be taken 
into account), transitions at smaller separations may be 
possible with appropriate materials. 
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Appendix A: Details of Computational Methods 

1. Scattering from periodic arrays 

There are numerous computational and semi- 
analytical methods from classical electromagnetic 
scattering that can be adapted to Casimir calcula- 
tions [52]. For example, there are computational 
techniques based on generic grids/meshes, e.g. finite- 
difference methods [32] [33] [53j [54] . There are also 
Casimir methods [28 , 29] based on classical cross-section 
methods (rigorous coupled- wave analysis [55], also called 
eigenmode expansion [56]), which divide the geometry 
into slices with constant cross-sections and match modal 
expansions between slices. Alternatively, there are 
spectral integral-equation methods tracing their roots 
to classic Mie scattering and related problems [57] : one 
divides the geometry into high-symmetry objects like 
spheres and cylinders, computing the scattering matrix 
for each object in a specialized basis (e.g. spherical 
waves), and then combining the matrices from different 
objects to match the boundary conditions. These 
methods have been adapted to Casimir problems for 
geometries consisting of a finite number of isolated 
objects [T6] [34H36] and corrugated surfaces [58]. Here, 
we adapt similar methods to periodic arrays of isolated 
objects (cylinders) by exploiting classical lattice-sum 
scattering methods [15] . For completeness, and because 
our formulation is directly in imaginary frequency 
and our normalization conventions differ from other 
applications p~4j [15], we outline the derivation in this 
section. 

The scattering matrix Tq from a single circular cylinder 
can be computed analytically in the basis of cylindrical 
multipoles, using the basis functions [T6] : 



E) p fe : re n g/out (^,x) = v < g r"(^x- Xj )z 



,reg/out 



(Al) 



where the linear differential operators L p (jp = M, N are 
the transverse electric and transverse magnetic polariza- 
tions, respectively) are: 



L M 



1 



i 



V x 
V x Vx 



and the cylindrical wave functions are: 



lOUt 

D k z ,n 



I n and K n are the modified cylindrical Bessel func- 
tions, and Xj gives the coordinate origin. The calculation 
of To is straightforward and we omit it here. We note 
that although the present method is simplified by the 
semi-analytic calculation of the scattering from circular 
cylinders, methods exist that can compute the scatter- 
ing matrices of non-circular cylinders of arbitrary cross- 
section [31]; such a hybrid method allows for efficient 
scattering computations from general two-dimensional 
arrays of compact objects. 

We now derive the scattering matrix T for a periodic 
array of cylindrical scatters in a uniform medium (the 
fluid interface will be added later). T is defined to be 
the scattering ma trix from each cylinder j in the basis 
given by Eq. (Al) with origin = jay, in the presence 



of the entire array of cylinders. This array has period a 
in the ^-direction and is translation-invariant in z, and 
only a single layer in the x-direction. In other words, 
for an incident field |E) m , the scattered field emitted by 
currents on cylinder j is: 



|-g(j)\scatt _ e ik y3 a Y\Ei) 



(A2) 



Here we have assumed that the incident field |E) m is 
a planewave with transverse wavevector components k y 
and k z , and have normalized it with respect to the origin. 
Then using the tran slation matrices [16 for the basis 
functions of Eq. (Al), the scattered field from object j, 



converted to the basis at the origin (j = 0), is given by: 

n" ,n' 

(A3) 

where the elements of the lattice sum S are given by: 

Sg = K m (|x| VkfT^e) e~ im9 > (A4) 

and 0j = sign^')^; and this transformation is diagonal 
in polarization p. Using linearity of the scattering pro- 
cess to sum over all incident cylindrical waves, the total 
incoming field at the origin xq for a planewave |E) m is ob- 
tained by summing the incident fields and the scattered 
fields from all cells j ^ 0: 
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IE' 



in, total 



1 + ( £ e ik y ja $ (j > ] AT 
(1 + SAT) |E) in 



|E) i: 



The matrix § is formally defined as the infinite sum of all 
S( J ) (in imaginary frequency this sum is well-defined due 
to the exponential decay of K n ), and A n ^ n > = (— l) n 5 n ^ n >. 
The field scattered from the cylinder at x = xo, in the 
presence of the array of cylinders, is then: 



| E y=o))Bcatt = T (l + SAT)|E) in 
= T|E) in 

where the last equality follows by the definition of T. 
Because this equation holds for arbitrary incident field 
vectors |E) m , the equation holds for the operators as well, 
allowing us to solve for T: 



T = [1 - T SA] 1 T 



(A5) 



Although both S and A are diagonal in polarization p, in 
general To is not and di ffere nt polarizations will couple to 
the periodicity via Eq. (A5). The scattering matrix from 



the entire array is a sum of T over all unit cells. This sum 
is more naturally expressed in a planewave basis, using 
the wave conversion matrices D for cylindrical waves [16] 
(here we have absorbed the normalization factors C p , de- 
fined in Appendix B. of Ref. \W\ into the definition of D), 
from which we find the planewave scattering matrix: 



In an anisotropic EM A, we assume that the i-th object 
is described by a homogeneous (position-independent) 
permittivity tensor 



e = 



e yz 



(A7) 



where the frequency dependence is left implicit and we 
used the fact that the permittivity tensor is symmetric 
on its spatial indices [59] . The relatively symmetric char- 
acter o f ou r setup leads to a very useful simplification 
in Eq. (A7) that becomes apparent in a coordinate sys- 



tem with its x-axis perpendicular to the slabs interfaces: 
it is clear that regardless of the in-plane orientation of 
the wires there is no mixing with the x-component, and 
so 



0. 



(A8) 



In addition, for our simplified EMA we make the further 



approximation that e x 



-yy 



s±. This symmetry im- 



plies that the e tensor is diagonal in the x, y, z coordinate 
system of Fig. [I] (top) and is of the form: 



e± 
e± 

£\\ 



(A9) 



Consider a planewave of transverse wavevector = 
(k y ,k z ) impinging on the slab, and let (j) be the angle 
makes with the z-axis. For our purpose, the scattering 
problem is best solved in the coordinate system (x f , y' \z') 
of the plane wave, where x' = x and z' is parallel to k^. 
In this basis, the e tensor becomes: 



r TB- 



K 



(A6) 



Here k is the reduced Bloch vector, and the matrix Rk 
couples all vectors k + miry for all integer m. 

The incorporation of multiple interfaces in the in- 
direction is achieved via a standard transfer matrix ap- 
proach, e.g., Ref. [15] For the case of a uniform medium 
£3 outside of the plates, the Fresnel formula are combined 
in a straightforward manner with the matrix T above to 
finally give the full scattering matrix R for the objects 
of Fig. g} 

For numerical computations, ten cylindrical multipole 
moments in the single-cylinder s catt ering matrix, 10/£ 
terms in the lattice sum of Eq. (A5), and a simple ex- 
ponential extrapolation from only 7x7 Brillouin zones 



for Eq. (A6) were sufficient for < 1% error. 



2. Effective- Medium Approximation 

In this section we outline the method used to predict 
the scattering matrices R from a specified e . The method 
is given in full detail in Ref. [TT1 



£ = 



e±e 
{s ir 







( £ \\ ' 
£||C 2 





- s_\_)sc 



where c = cos (j) and s = sin (j). 

We can now proceed to determine the scattering ma- 
trices of planewaves for the isotropic-anisotropic flat in- 
terface. After a considerably long algebraic calculation 
(for details see Ref. [T7J but note the change in coordinate 
systems between our Fig. [l] and their Fig. 3, and that for 
us c = cos(/> and the speed of light is set to 1), it is pos- 
sible to show that the four reflected/incident amplitude 
ratios are given by 



r TE ' TE (^,k ± ) 
r™' TE (^,k ± ) 
r TE '™(^,k ± ) 



detl 



detM 
detM 2 



detM 
det M 3 



JM.TM 



detM 
det M 4 



detl 



10 



where 



-1 







a( 2 ) 1 









-/3( 2 ) 





q in /u 


1 


1 





-1 


(-\ \ 




i 







Q ( J ) 1 


q in /uj 







_^(2) 





<2in/^ 


1 


1 





-1 






-1 


1 




«(2) 1 




<7in/^ 




_^(2) 








1 


1 








(A \ 












«(2) 1 










-/?(2) 






1 


1 


i 


-1 




f9^ 


-l 







«(2) ] 







-/?« 


-/3(2) 





q in /u 


1 


1 





1 


7 (i) 


7 (2) 



\ G y G z) 



( G y' G z 



\ G y G z) 



\G' y G' z ) {G y G z \KlG' y G' z 



FIG. 8: Feynman rules for the diagrammatic expansion of the 
Casimir energy, following Ref . 1451 In our notation, planewave 
states are indexed by their reciprocal lattice vector = 
(G y ,G z ); frequency, polarization, and the conserved reduced 
Bloch vector are suppressed as our focus is on non-specular 
(i.e., not conserved) reflections. Here U is the free-space 
translation matrix, and 1U the scattering matrix for object i 
(= 1, 2 for our case). Primes on the polarizations (e.g., y') 
denote distinct wavevectors. 



Finally, the elements Lij , are defined in terms of the per- 
mittivity tensor as 



U0 2 £±_ 



1 



^23 — 1 

L 31 = -L 42 (e± - e\\)sc 

UJ 2 

-\\ c 



^32 



• e±c 2 + e\\s 2 



L41 = -s\\(? — e±s 2 



A careful analysis shows that the equations Eq. (A10) 
become singular at (j) = nf for integer n. In this case, 
we take the limit <j) — >• n\ numerically. 



Appendix B: Suppression of non-specular scattering 



Here q in = ^u 2 - k^_; 
above are defined by 



\J uS 1 — ; the coefficients in the matrices 



(q^) 2 -u 2 A 
u 2 d 

L31 L% 



7 



Cp) 



L41 



L42 ( p ) 
7 (p) q(p) 



(A10) 



where p = 1 , 2 and also 

qM = q 
with 



{p) = n(±) = uj^-^A^- B ± y/(A - B) 2 + 4C (All) 
V 2 



A = L14L41 
^ = L23I/32 

C = ^14^42^23^31 



In this section we give a proof that non-specular 
scatterings (i.e., those not conserving the perpendicular 
wavevector k^ + Gj_) are greatly suppressed for crossed 
slabs relative to aligned slabs. This explains the wide 
range of validity for the specular-reflection approxima- 
tion for crossed slabs. The proof is a straightforward 
application of a recently-developed diagrammatic expan- 
sion for the Casimir force [45]. In this framework, the 
logdet expression of Eq. ^ is re-expressed as Tr log and 
expanded in a power series. Each term in the power se- 
ries can be computed via a set of Feynman rules, shown 
in Fig. [8] For sufficiently large separations, this series is 
rapidly convergent. 

Each scattering event is represented as the process 
\G X ,G Z ) — » \G' X ,G' Z ), where G X1 G Z represent the per- 
pendicular wavevector of the planewave. We suppress 
both the frequency and reduced Bloch vector (which 
are both conserved in the scattering process), and the 
polarization p (which is not relevant to our analysis of 
non-specular reflection). Each object is represented as 
a pair indicating the object alignment (Y, Z) and object 
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Wires Aligned Wires Crossed 




FIG. 9: Lowest-order diagrams representing the contribution 
of non-specular scattering events to the Casimir energy, con- 
structed using the rules of Fig. [8] Because we are measuring 
the energy of interaction between two objects, successive in- 
teraction vertices must involve different objects. On the left, 
when both objects 1 and 2 are directed along the z-axis, a 
non-specular event can occur in a 2-vertex diagram. How- 
ever, when the wires are crossed, the lowest-order allowed 
non-specular process involves a 4-vertex diagram. 



index (1,2). The key point is that, due to translation- 
invariance, the scattering matrix for a Z object is pro- 
portional to S (G z — G' z ), and for a Y object is propor- 
tional to S (G y — G'y). Further, the free-space propaga- 
tor is diagonal in the planewave basis and thus conserves 
both G y and G z . Successive interaction vertices in a dia- 
gram must involve distinct objects as external insertions; 
therefore, for the crossed configuration, a planewave must 
scatter twice as many times to return to its original state 
in the lowest-order diagram relative to the aligned con- 
figuration. Example lowest-order diagrams are shown 
in Fig. [9] We find that for aligned (e.g., ZZ) slabs, the 
lowest order diagram involving non-specular reflection in- 
volves two vertices, while for crossed slabs it involves four 
vertices. Since each propagator represents an exponen- 
tial attenuation, and each scattering event multiplication 
by a number of magnitude less than unity, this implies 
that non-specular interactions are greatly suppressed for 
crossed slabs. 
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